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Abstract. Field-theoretic construction of functional representations of 
solutions of stochastic differential equations and master equations is re- 
viewed. A generic expression for the generating function of Green func- 
tions of stochastic systems is put forward. Relation of ambiguities in 
stochastic differential equations and in the functional representations is 
discussed. Ordinary differential equations for expectation values and cor- 
relation functions are inferred with the aid of a variational approach. 
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^ • 1 Introduction 



Fluctuation effects in physics, chemistry, biology, operations research etc. are 
often described by the Langevin equation (stochastic differential equation, SDE) 

^ = -A> + C/(^)-(-/6M, (1) 

where / is a Gaussian random field with zero mean and the white-in-time cor- 
relation function 

{f{t,^)f{t',K'))^S{t-t')D{^-^'). (2) 

The Langevin equation with white noise is mathematically ill-defined. A consis- 
tent description may be obtained by using, e.g., a S sequence with finite corre- 
lation times 

{f{t, ^)f{t', x')) = D{t, x; t', x') ~^ S{t - t')D{^, x') , 

in the construction of the expectation values of functions of (p and passing to the 
white-noise limit in the end. This is often physically natural and leads to solution, 
which corresponds to the Stratonovich interpretation of the SDE. However, the 
point of introducing of the SDE with white noise is to avoid dealing with this 
limit explicitly. Technically, the Stratonovich interpretation is complicated and 
the SDE is most often used in the Ito interpretation in mathematical treatments. 

Instead of using the mathematically problematic, although physically trans- 
parent, Langevin equation the stochastic problem ([T]), ^ may be equivalently 
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stated in terms of the Fokker-Planck equation (FPE) , which is an equation for 
both the conditional probabihty density p {ip,t\ipo^to) and the probabihty den- 
sity p {if, t) of the variable ip. The simple way to demonstrate this equivalence 
uses rules of Ito calculus [2], however, and I am not going to dwell on this is- 
sue here. Other seemingly simple methods use - at least implicitly - functional 
integrals, which are still mathematically ill-defined, let alone the Gaussian inte- 
gral. Therefore, only the correspondence between the quantities specifying the 
stochastic problem in both approaches will be quoted here. The main advan- 
tage of the Fokker-Planck equation is that the equation itself is completely well- 
defined partial differential (or functional-differential for field variables) equation. 
The ambiguity of the Langevin problem shows in that the FPE is different for 
different interpretations of the SDE. 

For simplicity of notation, consider zero-dimensional field theory. The Fokker- 
Planck equation for the conditional probability density p {kp, i|<^Oi ^o) in the case 
of the Ito equation is 

d d 

—p{ip,t\ipo,to) = -—{[-Kip + U{ip)]p{ip,t\ipo,to)} 

+ 1-^ H^)Dbiip)piip,t\ipo,to)] . (3) 

If the SDE ([T]) is interpreted in the Stratonovich sense, the FPE is 
d d 

—p (if, t\Lpf), tQ) = ~— {[-Kip + U{ip)] p {(f, t\(po, to)} 

+ {^(^)^ [Db{^)p{p,t\^oM)]^ ■ (4) 

In what follows I will use the fact that the conditional probability density 
p{ip,t\ipQ,to) is the fundamental solution of the FPE ^ or i.e. 

p{Lp,ta\ipo,tQ) ^6{ip-ipo) . 

Contractions are not quite obvious, when the random variable has several com- 
ponents. For instance, the Fokker-Planck equation in the Ito form becomes 

d d 

—p{ip,t\ipo,to) = {[-Kijpj + Ui{(p)]p{ip,t\'po,to)} 

1 92 

+ na a [bik{ip)Dklbjl{(p)p {ip,t\ipo,to)] 

2 0(piO(pj 

for a multi-component variable ipi. 

The Fokker-Planck equation is similar to Schrodinger equation. Using this 
analogy, the solution of the FPE as well as calculation of expectation values 
may be represented in a way analogous to quantum mechanics [7] . Construction 
with the FPE as the starting point gives rise to the famous Martin-Siggia-Rose 
solution of the SDE [5], but avoids ambiguities inherent in the SDE. 
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Markov processes described in terms of the Fokker-Planck equation have 
continuous sample paths, i.e. = (p(t) is a continuous function of time. Not 
all interesting stochastic processes belong to this category. A wide class of such 
processes describe changes in occupation numbers (e.g. individuals of some pop- 
ulation, molecules in chemical reaction) which cannot be naturally described by 
continuous paths. This kind of processes are described by master equations - a 
special case of (differential) Kolmogorov equations ^ . 

The generic form of a master equation written for the conditional probability 
density p (ip, t\ipQ, to) of a Markov process is 

d f 

—p{ip,t\ipo,to)^ dx[W{ip\x,t)p{x,t\ipo,to) -W{x\'P,t)p{ip,t\ipo,to)] , 

where W^(iy9|x,i) is the transition probability per unit time, whose formal defini- 
tion from the differential Kolmogorov equation is (for all e > 0) 

T^(^|x,0=Jjm , 

uniformly in ip, x S'Hd t for all — xl > £• 

I shall be using the master equation for discrete variables (occupation num- 
bers). In this case the discontinuous character of the paths of the jump processes 
described by the master equation is especially transparent. The transition proba- 
bilities are usually simple functions of the occupation number n. As an example, 
consider the generic master equation for the stochastic Verhulst model 

^E^hA = [^(^ + 1) + + lf]P{t, n + l) + X{n - l)P{t, n-1) 
at 

- {Pn + \n + 7n2) P(t, n) , (5) 

where /? is the death rate, A the birth rate and 7 the damping coefficient necessary 
to bring about a saturation for the population. The choice X = (3 — Q leads to 
the master equation for the annihilation reaction A + A ^ A. The set of master 
equations may also be cast in the form of an evolution equation of the type 
of Schrodinger equation in a Fock space of many-particle quantum mechanics. 
This representation is due to Doi [T]. Contrary to the widely known Martin- 
Siggia-Rose approach, the method of Doi has been gaining due attention only 
recently. 

The two quantum-mechanical treatments of stochastic problems have rather 
different appearances in the literature. The first aim of this paper is to present 
both cases in a unified fashion in terms of generating functions of Green func- 
tions of many-body quantum mechanics. Second, ambiguities in the functional 
representation of the solution of the stochastic problems are expressed in terms 
of the thoroughly analyzed functional form of the quantum field theory. Third, 
the functional representation is used to calculate fluctuation corrections to rate 
equations. 
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This paper is organized as follows. In Section [5] the operator representation 
of stochastic problems and their connection to Green functions is reviewed. Sec- 
tion |3] is devoted to discussion of functional representations of Green functions. 
Variational method in the functional approach is discussed in Section |4] and 
conclusions presented in Section [S) 

2 Quantum-Mechanical Representation of Stochastic 
Problems 

2.1 Green Functions for the Fokker-Planck Equation 

Consider, for definiteness, the Fokker-Planck equation ([3]) corresponding to the 
Ito interpretation of the Langevin equation ([T]). Introduce - in analogy with 
Dirac's notation in quantum mechanics - the state vector | pt ) according to the 
following representation of the PDF 

which is the solution of the FPE (jS]) with the initial condition p{ip, 0) — po{(p)- 
To construct the evolution operator for the state vector, introduce momentum 
and coordinate operators in the manner of quantum mechanics by relations 

d 

In these terms, the FPE for the PDF gives rise to the evolution equation for the 
state vector in the form 

— |pt) =L\pt), 

where the Liouville operator for the FPE corresponding to the Ito interpretation 
of the SDE assumes, according to the form 

L = Tr[-K^ + U{ip)] + ^n^b{ip)Db{0). (6) 

Note that, contrary to quantum mechanics, there is no ordering ambiguity in 
the construction of the Liouville operator here. 

In this notation, the conditional PDF may be expressed as the matrix element 

p{ip,t\ipo,to) = l^ip 

Introduce time-dependent operators ip{t) in the Heisenberg picture of imaginary 
time quantum mechanics (i.e. Euclidean quantum mechanics): 



(7) 



(8) 
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and define the time-ordered product (chronological product, T product) of time- 
dependent operators 



T 



Ai{ti) ■ ■ ■ An{tn) ^ ^ P e [h . . .tn) Ai{h) ■ ■ ■ A^{tn) , (9) 
P(l,...,n) 



where 

e{h... i„) = e{h~ t2) 9{h-h)---e (t„_i - i„) . 

In definition ^ the sum is taken over all permutations of the labels of the 
time arguments {ti}"^^ and the operators in each term are put in the order of 
growing time arguments from the right to the left. Thus, under the T-product 
sign operators commute. It should be noted that the definition of the time- 
ordered product should be amended for coinciding time arguments. 

Introduce then the n-point Green function as the quantum-mechanical ex- 
pectation value of the T-product of n operators ([5]) 

G„(ti,i2,...i„) =Tr {poT [ipHih)ipHit2) ■ ■ ■ 0H{tn)]} (10) 

determined by the trace Tr and the density operator 

Po^ dy:>\po){(p\. (11) 



Choosing, for definiteness, the time sequence ii > ^2 > ^3 > • • • > in-i > tn > io 
it is readily seen by direct substitution of relations ([7]), ^ and pT|) in pUjl 
with the aid of the normalization conditions of the PDF and insertions of the 
resolution of the identity J dip \ip){ip\ = 1 that 



d(pi ...J dipn (/?! • • • ipnP{'Pl,ti; ip2,t2] Vn.tn) = Gn{tlMi ■ --Ui) , (12) 

i.e. the Green function pUj) is equal to the moment function ([T^ . This relation 
connects the operator approach to evaluation of moments of the random process. 

2.2 Green Functions for the Master Equation 

The set of master equations for P{t, n) is reduced to a single equation by "second 
quantization" of Doi. Build a Fock space: operators a, d^ and basis vectors | n ): 

a|0)=0, d^|n) = |ri + l), [a, a+]==l, {n\m)~n\6nm- 

Master equations yield kinetic equation for the state vector 

oo 

in the form of a single evolution equation for the state vector (jl3p without any 
explicit dependence on the occupation number: 

^=L(a+,a)|P*), (14) 
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where the Liouville operator L{a^ , a) is constructed from the set of master 
equations according to rules: 

nP{t, n)\n) = d^aP{t, n)\n) , 
nP{t, n) \ n — I ) = aP{t, n)\n) , 
nP{t,n)\n+ 1) = a+a'^dP{t,n)\n) . 

For instance, the Liouville operator for the stochastic Verhulst model ([5]) is 

L{a+,a) = f3{I - a+)a + -f{I - a+)aa+a + X{a+ - I)a+a . (15) 

Equation (jl4p gives rise to the Heisenberg evolution of operators in analogy with 
([5]). To represent expectation values of occupation- number dependent quantities 
in the operator formalism use the projection vector (P |: 

oo oo 

^ — ' ni ^ — ' n! 

Consider the Green function of occupation- number operators fiH{t) — d'^{t)dH{t): 

Gm{tl,t2, . . . im) = Tr [PoT{nHiti)nH{t2)---nHit,n)]} , (17) 
with the density operator 

Po^\Po){P\. (18) 

From definitions it follows that the conditional probability density function for 
the master equation may be written as (the factorial in front of the matrix 
element is due to the unusual normalization of the basis states) 

P{n,t\no,to)^^ (n e^(*~*°) m) . (19) 



n 



Choosing, for definiteness, the time sequence > ^2 > ^3 > • • • > ^n-i > tn > to 
it is readily seen by direct substitution of relations ([T9| and (fT8)) in (fT7|) with the 
aid of the normalization conditions of the PDF and insertions of the resolution 

of the identity 7^ — I n ) ( n I = 1 that 



^ . . . ^ ni • • ■n„iF(ni,ii; 712,^2; • • • ;ri„i,t„i) = Gm{ti,t2, ■ ■ .1^) , (20) 



i.e. the Green function P7|) is equal to the moment function (|20p. This relation 
connects the operator approach to evaluation of moments of the random process 
described by a master equation. 
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3 Functional Representatation of Green Functions 

The generic representation for the generating function of Green functions intro- 
duced in the previous section may be written in the form 



G( J) Tr 



where the source term and the density operator are 



dt(pH{t)J{t) , po = d^\PQ){v\ 



Sj 



for the Fokker-Planck equation, or, for the master equation. 



(21) 



Sj 



dta+it)aH{t)J{t), /}o = |Po)(P| 



Perturbation theory is constructed in the Dirac (interaction) picture of quantum 
mechanics. To this end, the Liouville operator is decomposed to the free and 
interaction parts Lq and Lj: 

L = La + Li. 

In our cases, the convenient choices are Lq = —K-ifip for © and Lq ~ —/3a~^d for 
(fT5|) . where K > 0, > 0. Time-dependent operators in the interaction picture 
are defined according to 

The corresponding evohition operator may be expressed in terms of the T- 
product (here, t > t' > to): 



U{t - to,t' - to) = e-(*-*o)ioe(*-*')^e(*'-*'')^'' 



= e^°*°rexp 



Liiu) du 



(22) 



Note that the Dirac operators in the T-product do not carry any ^o-dependence 
in ([22]). In the interaction picture the T-product in (j2ip may be written as [9] 



5^°*«?7(io,i/)r 



U{t,,to)e'''<>'° , 



(23) 



where Sj — J^' Lj{t) dt and the time instants tf and ti are chosen such that tf > 
ti > ti > to for ah ti, I = 1,2, ... ,n. The evolution operator with the reversed 
order of time arguments may be cast into the form of the anti- chronologically 
ordered exponential (t > to) 



(24) 
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With the use of representations (|22l) and ([M)) the chronological product in the 
generating function (1^1]) may be expressed as the product of three chronological 
products: 



Te' 



T 



.Sj+Si 



Te 



Ito Lit) dt^ Loto 



(25) 



The three T-products in ([25)1 fuse - due to Wick's theorems - in a normal product 
(tt to the left of cp or a+ to the left of a) giving rise to the representation [5] 



G{J) = Tre-'^"*»poe 



X exp 



S. 



JW2) 



exp 



to 



2 6^ ~ 



bi) du - 



15^5 
1 ' 2 5(1)2 S(p 
tf 



_6 6_ 

54>i 5cj)2 



Li{4)2) du 



f>l=02=</' 



(26) 



where is a two-component shorthand for all the operators appearing in Li, i.e. 
either (p — (i^,7t) ot (j) = (d,a+). In (1261) . the auxiliary field variables 0i and (j)2 
correspond to functional arguments prescribed to the anti-chronological product 
and the chronological products in (|25p , respectively. Originally, in the functional 
representation of Wick's theorem for (1251) to each T-product a separate field 
variable is prescribed, but those corresponding to the two consecutive rightmost 
T-product factors may be replaced by a single variable (^2 ■ The propagators and 
contractions in (pS)) are standard: 



A{t,t') = T mkt' 



The functional Lj in (j26p (not an operator any more) is ambiguous, however, 
because the chronological products have not been defined at coinciding time ar- 
guments. The choice of the value of the T and T products at coinciding time 
arguments fixes the form of the interaction functional Lj. It should be noted 
that this choice affects the propagators as well. The normal form obtained by 
replacing operators by functions in the representation obeying Lj = N[Li] (see 
[S] for details) is the most natural here. For the T-product this boils down to 
choosing the temporal step function equal to zero at the origin. Thus, the practi- 
cal rule to resolve this ambiguity sounds similar to that in the Ito interpretation 
of the SDE [2], but it should be borne in mind that that the latter ambiguity 
has already been fixed by different means and there is no obligation to make the 
same choice here. Henceforth, the normal form of Lj is implied together with 
the corresponding choice: A{t,t) = 0, A{t,t) — 0. 

To cast the generating function in the functional form without operators, the 
expression 

TTe~^°*"poe^°*°N[...] (27) 
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should be calculated. Here, the ellipsis under the normal product sign denotes 
an arbitrary operator. This calculation is different in the two cases discussed. 
Moreover, in practical calculations for the master equation it is convenient to 
pull the coherent state exponential of the projection vector ((T6)) through all the 
operator expressions to the rightmost ket- vector of the state vector ([T5| . This 
leads to the shift of all creation operators by the unit operator: a'^ ^ d'^ + 1. 
I shall keep the free Liouville operator in the form chosen above and put all 
changes in the interaction operator: L — > Lq + L'j. After these manipulations the 
density operator in ((27)) for the master equation case assumes the form (strictly 
speaking, this is not a density operator any more) 

oo 

P,^\P^){0\ = Y^P{Q,n) [a+ + l]"|0)(Ol. 

Tl = 

It is an important property of the field theory of these stochastic problems 
that in representation ([26l) the terms corresponding to the chronological product 
(expressed in terms of the field ^2) and the anti-chronological product 
effectively factorize and the latter gives rise to the unity. This is different from 
the Green functions of quantum kinetic theory based on a similar expression. 
The factorization of the anti-chronological part in (051) is a consequence of three 
features. First, only contractions n{t, t') = ip{t)n{t') — N [ip{t)n{t')] and n{t, t') = 
d{t)a^{t') — N [a{t)a^{t')\ are non-trivial, all the rest vanish; second, the average 
(P7|) does not add anything to propagators and, third, in the interaction terms 
the two fields 0i, 4>2 are not mixed. 

Results of calculation of ([26]) will be quoted here for two important cases: 
First, the probability distribution is assumed to be homogenous, i.e. that there 
exists an equilibrium limit 



p{ip,t\ipo,to) Pei^p) , to ^-00, 



second, the limit to^ti — > —00, which allows to consider the problem on the whole 
time axis. 

With the use of the standard manipulations |3| in the calculation of (|27l) I 
obtain in the case of the Fokker-Planck problem the representation (with the 
choice i/ — ?> 00) 



G(J) = / dr]pe{ri) exp 



S(p Stt 



X exp 



[Li{t) + ip{t)J{t)] dt 



■K = 



(28) 



where (/?^ is the solution of the equation (dt + K)Lp,j{t) — with the initial 
condition ipn{0) = rj. Here, the choice ti = Q gives rise to the standard Cauchy 
problem, whereas in the limit ti ^ ~ 00 leads to the "scattering" problem on the 
whole time axis with the limit (p^ — > 0. In that case the normalization condition 
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oi Peiv) yields 

G( J) cxp 



—A— 

Sip Sn 



exp 



[Li{t) + ip{t)J{t)] dt 



1^ = 
7r = 



(29) 



which, in fact, is the standard functional representation of quantum field theory 
(QFT) 9 and allows to use the methods of the latter in practical calculations. 

In the case of the master equation problem the corresponding functionals for 
the generating function are 



G(J) 



ds 



2m 
X exp 



e''G(s) exp 



6a Sa+ 



where 



G{s 



L'j{t) + [{a+{t) + l)a(t)] J{t)] dt 

E°° r(n) 



i+=0 



(30) 



and fls is the solution of the equation {dt + f3)as(t) — with the initial condition 
0^(0) = s. The generating function for the problem on the whole time axis 
assumes the form 



G(J) = exp 



6a (5a+ 



X exp 



[L'j{t)+[{a+{t) + l)ait)] J{t)] dt 



a=0 
i+=0 



(31) 



Although these expressions are written for the zero-dimensional field theory, 
their generalization to fields in finite-dimensional space is straightforward. In 
that case K and (3 stand for positive-definite second-order differential operators. 

Functional-integral representations for generating functions may be obtained 
with the aid of the Gaussian integral representation 



exp 



S ^6 

6ip 6t: 



detA / / Va+Va exp 



-a+A ^a + a+— + a— 



6lp 



(32) 



of the functional differential operators (for brevity, only the case of the Fokker- 
Planck equation will be discussed here). It should be noted that the kernel A~^ 
on the right side of p2l) contains a differential operator. Therefore, the functional 
integral generated with the aid of (15^ must imply boundary conditions in the 
space of integration which single out the proper Green function of this differential 
operator as the propagator A on the left side 0. 

In addition to this, the choice of the value of the chronological product at 
coinciding time arguments affects both the form of the interaction functional 
/ Li(t)dt and the value of the propagator A at coinciding time arguments. As 
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discussed above, this effect is explicit in the interaction functional, although often 
overlooked. It is not clear, however, how this choice may be taken into account in 
the Gaussian integral (l32|) . The ambiguity in the definition of the chronological 
product bears resemblance with the Ito-Stratonovich ambiguity of the original 
Langevin equation, because its solution boils down to the choice of the value of 
the step function at the origin. It should be borne in mind, however, that I have 
got rid of the latter at the outset by starting from the Fokker-Planck equation. It 
is worth recalling that at any rate the perturbation expansion is independent of 
the choice of the amendment of the definition of the chronological product and 
this feature must be reproduced by any alternative calculation of the functional 
integral for generating functions. 

However, in the functional-integral representations emerging from (j28p . (j29p 
and (j32|) ambiguities inherent in the functional solution of the stochastic problem 
are explicit only in the interaction functional. Therefore, care has to be exercised 
to take into account subtleties discussed above, when the functional integral is 
used for numerical calculation of the Green functions - or for any calculation 
other than the perturbation expansion. The main point of the formulae ([25)) - pip 
is to give a unified, compact and unambiguous representation of the perturbative 
solution of the stochastic problems discussed. 



4 Variational method 



From the wealth of calculational methods of the QFT useful in stochastic prob- 
lems [lOj I would like to emphasize here the variational approach. The more 
or less standard steepest descent method has been applied to calculation of 
the functional integrals for generating functions, e.g., in order to evaluate the 
asymptotic behaviour of high orders of the perturbation expansion. Since the 
method is based on the knowledge of an explicit solution of a nonlinear differen- 
tial equation, useful results have been scarce. Apart form the zero-dimensional 
field theory, only few cases - based on the previously known instanton of the 
static tf^ model - have been analyzed for various models of critical dynamics 



The variational method may also be applied to construct fluctuation-amended 
rate equations and the like. To this end, the generating functional G is not suit- 
able because it does not give rise to solutions in terms of moments directly. Such 
solutions may be obtained with the aid of the Legendre transforms of In G with 
respect to sources. To obtain rate equations, it is convenient to construct the 
generating functional (only the master-equation problem will be discussed here) 
with sources linear in a"*" and a: 




oo 



5^P(0,n) [a+(0) + l] 
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Here, S'(a+,a) is the dynamic action determined by the differential operator 
giving rise to the propagator and the interaction functional Vf. 

S{a+,a) = Jdt{a+ [-~dt ~ (3] a + L'j{a+ ,a)} . 

The functional Legendre transform is defined via equations 

SlnG ^ SlnG 



r{a+,a) = \nG{J+,J) - a+J+ ~ aJ , 



SJ ' 5J^ 



where the (iterative) solution of the latter with the subsequent substitution in 
the expression for /^(a"*", a) is implied. Seek now a"*" and a as the extrema of the 
functional 

r(a+, a) = r(a+, a) + a+ J+ + aj . 

This approach has the advantages that the problem is formulated in terms of 
explicit functions of the first moments a"*" and a of the random fields and it 
allows to find possible nonperturbative solutions. 

Unfortunately, no closed form is known for the functional r{a~^ , a). The rules 
for the construction of the perturbative series are, however, well known [9]. The 
corresponding expression for the functional F(a^,a) starts with terms corre- 
sponding to the dynamic action and the initial condition, the rest consisting of 
an infinite sum of Feynman graphs corresponding to the terms of perturbation 
expansion: 

F{a+,a) ^ S{a+,a) + a+J+ + aJ + \nl ^ F(0, n) [a+(0) + l]"l + loops. 

ln=0 J 

For instance, for the reaction A + A with the initial Poisson distribution 
in the second order in the coupling constant A the effective action is 

r = - ^ dtjdx {a+dfa - Da+W'^a + \D [2a+ + [a+f] a^] 

+ nojdxa+{x,0) + ^ ■'HCH^ +1 ^Kin>^C +■••' 

where the wavy line denotes a and the wavy line with slash As a result 
from variations with respect to a and a"*", the fluctuation-amended rate equation 
follows in the form 



dta = DV^a - 2XDa^ + 4X^D^ du J dy A' 



(t-u,x-y)a2(w,y) + ..., (33) 



where A is the diffusion kernel {dt — DV^)A — 1. The analysis of the fluctuation 
term shows that (1551) is not well defined as it stands. The fiuctuation-amended 
rate equation should be written in terms of a renormalized field theory. In that 
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case the asymptotic behaviour of its solution may be analyzed with the aid of 
the field-theoretic renormalization group, which has been used for the analysis of 
the effect of density and drift fluctuations on the long-time decay already at the 
rate-equation level [6l3j . The analysis of the renormalized counterpart of p3|) is 
in progress. 

5 Conclusion 

Construction of functional representations of solutions of stochastic problems is 
reviewed. A unified generic representation of the generating function of Green 
functions of random fields brought about cither by the Fokker-Planck equa- 
tion (corresponding to a Langevin equation) or by the master equation is de- 
rived within the operator approach. Relation between the ambiguities due to the 
stochastic differential equation and the functional representation are discussed. 
Practically important cases of the Cauchy problem and the scattering problem 
are presented. Significance of the choice of the functional representation for nu- 
merical calculation of response and correlation functions is analyzed. Importance 
of the variational approach to evaluation of moments of random fields is empha- 
sized. The fiuctuation amended rate equation for the reaction A + A ^ is 
discussed. 
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